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Abstract 

With  the  assumption  of  constancy  of  the  elastic  constants 
in  Hooke ' s  Law  over  the  range  of  strains  involved,  but  no  limit- 
ation upon  magnitude  of  strains  or  displacements,  formulae  for 
strain  energy  are  derived  for  several  cases  of  interest.   Exact 
buckling  criteria  are  derived  from  the  effect  of  small  pertur- 
bations, from  displacements  in  an  unbuckled  state,  upon  the 
total  energy. 

This  method  is  applied  to  three  problems:   A  column  of 
rectangular  section,  compressed  so  as  to  keep  its  ends  flat  and 
parallel  to  their  ori-<lnal  position;  a  hollow  circular  cylinder, 
subjected  to  uniform  pressure  on  its  outer  surface;  and  a  hollcv. 
sphere,  subjected  to  uniform  external  pressure.   In  all  cases, 
for  thin  sections,  results  a~ree  with  those  derived  by  other 
means.   It  is  shown  that  vjith  the  given  assumption,  a  buckling 
or  critical  load  always  exists,  no  matter  what  the  geometric 
proportions .   As  an  interesting  sidelight,  it  is  also  shown  that 
a  finite  pressure  may  expand  the  cylinder  or  sphere  to  infinity. 
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1,   Introduction. 

Assuming  linearity  between  stress  and  strain  when  referred 
to  principal  directions,  an  exact,  non-linear,  formula  for  the 
strain  energy  may  be  derived  in  terms  of  displacements.   Althou^.i 
not  shown  in  this  paper,  the  same  i'ormula  results  if  Prof. 
Friedrichs's  form  of  Hooke ' s  Law  for  large  displacements  is 
utilized.   The  work  done  a^^ainst  external  forces  is  also  evalua- 
ted and  added  to  >7ive  the  total  potential  energy  in  the  system. 
Buckling  occurs  when  a  small  perturbation  of  displacements  cor- 
responding to  unbuckled  conditions  results  in  a  decrease  of  the 
total  energy.   Since  the  perturbation  may  be  taken  arbitrarily 
small,  higher  powers  in  the  displacements  corresponding  to  it 
and  in  derivatives  of  such  displacements  can  be  dropped  without 
error  leaving  only  quadratic  terms  to  be  considered,  linear  terms 
vanishing  since  the  energy  must  be  an  extremum  for  the  unbuckled 
state . 

The  change  in  energy  due  to  the  perturbation  is  now  maxi- 
mized by  application  of  classical  Calculus  of  Variations  rules 
to  yield  linear  differential  equations  determining  optimal  rela- 
tions for  the  perturbation  throughout  the  interior  of  the  body. 
These  are  solved,  subject  to  boundary  conditions,  in  terms  of 
arbitrary  constants  and  the  displacements  obtained  are  substitu- 
ted back  into  the  expression  for  change  in  total  energy.   The 
limits  of  positive-dof initeness  of  the  energy  x^;ith  respect  to 
the  integration  constants  give  the  criteria  for  buckling. 
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2 .   The  Strain  Energ,y  in  Terms  of  Principal  Straj.ns . 

If  e-,  ,  Bp,    e^   are  the  unit  strains  in  the  principal  di- 
rections at  any  point  in  a  strained  body  and  p,  ,  Pp,  p^  are  the 
corresponding  stresses  oi  forces  per  unit  area,  we  assume  that 
the  latter  are  linear  functions  of  the  forrner  so  that,  for  iso- 
tropic material,  we  have 

p^  =  \(e-|  +  62  +  e^)  +  2ij,6^  ;    i  =  1,2,3 

where   X   and   \i   ^re   elastic    constants    (Lame).      If  w    is    the   strain 
energy,   we   mast   also   have 

?)w 
Pi   =   BeT 

so  that  the  energy  must  be  given  by 

(1)  2w  =  \(e^  +  62  +  e^)^  +  2|i(e^  +  e^   +   e^) 

If  we   assume    it   to  be    zero  for   zero    strain. 
If  p^   =   0,    then 

®3   "    ■  X?2ll   ^®1   ^   ®2^ 
and   the   above    reduces    to 

(2)  2w  =  ^^^^-    (e^   +   e2)^  +   ix(e^    -   e^)^ 
If   e      -   Oj    then  we  pet   the    similar   form 

(3)  2w  =    (X  +   ^L)(e^    +   e^)^   +   ix(e^    -   62)^ 
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3 .   Principal  Strains  in  Temns  of  Displacement , 

Consider  a  two-dimensional  situation  In  Cartesian  coordi- 
nates.  A  point  with  coordinates  x,y  moves  to  position  x+u,y+v 
when  the  body  Is  strained.   An  adjacent  point  an  Incremental  dis- 
tance dS  away  at  angle  a  to  t  he  x  axis  has  Initial  coordinates 
X  +  dS  cos  a,y  +  dS  sin  a  and,  after  strain^ is  dS^  away  at  angle 
p  with  coordinates  x+u+dS,  cos  p,y+v+dS,  sin  p.   Assume  these 
two  po^'nts  to  deterralne  one  principal  direction  of  strain  (see 
Pig.  1).   The  second  principal  direction  is  at  right  angles  to 
the  first.   A  point  on   it  an  incremental  distance  dS  from  our 
first  point  has  initial  coordinates  x-d3  sin  a,  y+dS  cos  a.   In 
the  stialned  state,  the  principal  directions  remain  perpendicu- 
lar, but  strains  may  differ  so  that  the  coordinates  of  the  third 
point  become  x+u-dSp  sin  p,  y-f-v+dSp  cos  p. 

Since  the  points  were  taken  very  close   to  each  other, 
continuity  of  displaceTient  gives  the  relations: 

dS,  cos  P  =  (1  +  u  )dS  cos  a  +  u  dS  sin  a 
1     ^    '     X'  y 

dS,  sin  P  =  v  dS  cos  a  i-  (1  +  v  )d3  sin  a 

1       X  y 

dSp  sin  p  =  (1  +  u  )dS  sin  a  -  u  dS  cos  a 
c.  X  y 

dSp  cos  p  =  -v  dS  sin  a  +  ( 1  +  v  )dS  cos  a 

"^  "^  y 

Eliminating   a   and   p    from   these   I4.  equations   gives    the    following 
2  equations : 

dS^        dS^ 

T^T—    •    T-g—  =(l    +    u)(l«-v)-uv 
d:3  dS  X  y'  y   X 
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so  that 

and 

(5)  (1  _  _2^2  ^  (    _  g  ^2  ^  (    -  V  )^  +  (u  +  V  )^ 
•^'      dS    dS  '     ^1    2       X    y       y    x 

In  polar  coordinates,  a  similar  procedure  can  be  followed 
(see  Pig.  2)  to  give  the  continuity  relations: 

dS,    cos    P  =    (1   +   r,  )dS  cos    a  -  i-  r^dS   sin   a 
±  p  P     y 

dS,    sin   P  -, 

— ~i-?~~  ^    -  ^p^""  cos    a  +  -   (1  +  4g)dS   sin   a 

dSp   sin  p  =    (1   +   r    )dS   sin  a  +  —  r^dS   cos    a 
d.  p  P      ^ 

dSp   cos    P  -J 

— =— 7 =   4  dS    sin   a  +  -4    (1   +  c|A)d3   cos    a 

p   +   r  ^p  P  Q 

and   the    strain  equations 

(6)  (2+9^  +  92)2  -    ^-^-^-^f-^^Q^-^  4o)  +  rp]2  +  [(p+r)4p-i   r^J^ 

(7)  (e^-e^)^^    [(f-4e-^f  4^)    "    r/   .    [(p.r)4p   +  ^  r^J^ 


h •   The  Rectangular  Column . 

We  take  (Fig.  3)  the  origin  at  the  mid-point  of  one  edge 
with  X  in  the  direction  of  the  load  and  y  at  right  angles  and 
assume  zero  stress  in  the  third  direction.   We  restrict  any 
buckling  so  that  the  ends  remain  flat  and  parallel  to  their 
original  position.   For-  uniform  compression,  since  there  are  n-^ 
forces  perpendicular  to  the  column  axis, 
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so  that 


e,  =  -e 


u  =  -  ex 


Xe 

If  this  is  a  stable  condition,  any  slight  change  in  displace- 
ments u  and  V  must  result  in  an  increase  in  strain  energy.   Con- 
versely, if  there  exist  small  changes  in  u  and  v  which  reduce 
the  strain  energy,  the  condition  of  uniform  compression  is  un- 
stable.  We  are  thus  led  to  assume 

u  =  -ex  +  P(x,y) 

and  investigate  the  ch;'-.nge  in  the  total  strain  energy  due  to 
small  P  and  G.   Since  we  are  considering  no  change  in  column 
length  between  the  uniform  compression  and  the  buckled  states, 
there  is  no  work  done  against  the  excorral  load  that  needs  to 
be  considered.   We  thus  find 
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We  also  note  that  our  boundary  conditions  require  F  and  G  to 
vanish  at  x  =  y  =  0  and  that  F  is  also  zero  for  x  =  0  or  x  =  L . 
Substituting  the  above  into  equations  ([|.)  and  (5)  and  the  results 
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into    (2),    we   obtain 
2w 


a:  +   2tj,    W  2(Xi-:x;        X        y'  '    y        x' 


*^tF^-0y-<f^e)2.,(Py.G^)' 


or 


2w        liX+2iil  ^2   _   2(3Xi2ix)    ^^     +    (p    .  g   )2   +    (p    +  G   )2 
ji'^X  +  p,  X+ti  X  X       y  y       x' 

■    (3^+2m-)    /p    .  ^    ^2  (3X+2M,)e  ,p        „    .2 

"  X+2p..         ^^x        y^       ~   i4-(A+ix;  -  (X+^^tje    ^^y  "     x' 

retaining   terms   only    to   second   order   in   the  derivatives   of  P  and 
G.      The    total   strain  energy  per  unit  width  of   column  is 
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W  =   j  j       w   dxdy 
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Since 


,L 


F     dx  =   F 
'  X 


=    0 


0 


we  have 


where 


2AW 


[I  K    +    [I  jj. 


A^  =  Q  dxdy 

h  1  ^, 


and 


Q  =    (P    -  G   )2+  (P    +G    )^+^4~%^-  (P    +G    )2  . -_l3Vf  2u|  e^      ,  ; 

X       y'         ^    y        x'  X    i-  d[i        X        y'        k(  A+ti^ -( A+2^e   y     x' 

For  e   <  0,    Q  >  0   so    that   /\w  >   0  and   the   uniform  stress   condi- 
tion  (tension,    in  this    case)    is    stable   as   expected.      Let   us    take 


-       ^    -     .'H    in    + 


i      / 


•■i  '  '^ 


:-'-      -i- 


»!i:i   r:.   'tobfo    ' -'ooj;; 


.'..■■ic 


ft    r..r'->''j' 


t«"    < 


-   WA 


I     ^r      , 


iLi.  —li 


••^'-f  «■      !' 


-^■•\    7--J:t    r-    0    <    ^v    .0 


8 
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M  =  if -t -fix 


and   change  variables    to 


^         L 


then 

2 

^  Q  =    (l  +  M)(P^  +  Gf)  +  2(M-1)F'  Cr    +  (1-K)(pf.+  G^)  +  2(1+K)F   G^ 

Since  P  =  0  at  Q  =  0  or  n,  we  may  express  It  in  the  half -range 

Fourier  series: 

P  =  T~  P  sin  nQ 
> —  n 

where  the  P  are  functions  of  >? .   We  also  note  that  v  =  0  at 
n  (  X 

0  =  0  or  Tt  so  that  G  may  be  developed  in  the  half-range  cosine 

series : 

G  =  G  -i    y  G     cos  nQ 
o   ^=^—  n 

If  we  use  primes  to  denote  differentiation  vjith  respect  to  >i  , 
we  thus  get 

Pq  =  XI  "F^ri  ^°^    ^® 

F     ==  XI  P^  sin  nQ 
G^  = 


0  =  -  H  "^^n  ^^^  ^'^ 


G„  =  G'  +  y"  G'  cos  nQ 
V)    0  ^ —  n 

Because  of  the  orthogonality  properties  of  the  trigonometric 
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terms,   we  have 


2 

Q  dQ  =  I 

0  -'  0 


Qdx  =  |  QdQ  =  £  (R^   .HR^) 


where 

R     =   2(1   +  M)G'^ 
o  o 

R^  =    (l+M)(n^P^+G^2)  +  2n(M-l)P^G;^+  (1-K)(f;^2+  n^G^)-2n(l+K)P.;^G^ 
Since   R      is   positive,    we  must   have   G'    =   0   for  minimiom  energy. 


which  makes 


G,       -       -      1.      G^        ly.O 


For   any  n,    /\W   is    an  extremum  if  we   satisfy   the    conditions 


or  when 


and 


9R  $R      , 

— ^  =    (—2:) 
n  n 


S>G  ^dG'^ 

n  n 


n(K  +  M)G'    =    (1   -  K)P"    -  (1   +  M)n^P 

■  n  n  n 

■n(K   +   M)P'    =    (1   +   M)G"    -  (1    -   K)n^G 

■  n  n  n 


from  which 


P  ^^-   2n^P"   +   n^P     =   0 
n  11  n 


or     P     =  A      sinh  nri  +   B  nri  cosh  ii);    +  0      cosh  ny)  +  D  nri  sinh  nr) 
n         n  '  n    '  '  n  '  n    '  ' 

and,    substituting   this   back   into   the  equations: 
^n  =    ^-^n^^^M^  B^Jcosh  nr,  -  B^nr)   sinh  nn 

■"    ^-^n^  ^SiP-  D,Jsinh  n>;    -  D^nr^  cosh  nr^ 
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Using   the    result    of  Appendix  A,    this  makes 


2 


H 
2 


R     dyl  =    [  (l-K)P  F'  +  (1+M)G  G'  +  n(M-K-2)P   G   ] 
ni  'nn  'nn  'nn 


H 
2 


which  upon  substituting  and   simplifying  yields 


H 

2 

D  ^v.   -    a2   sinh  nH      »    ^ 

I  — 77  1     „  R      a>i    =   A      g A   B 

L|.nH  J    H     n  (          n       nH              n  n 
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nH 


2(M+K)' 


(2+M-K)    sinh  nH  ,       S] 

- — --  -    f   ^^ cosh  nH 

.  .       K) 
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M^K  nH 

(1+M)(1-K)(2+L-K)    sinh  nH      (1+M)(1-K) 


(1+MK) 


nH 
2  nH 


cosh  r''" 
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(2+M-K)    sinh  nH 
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cosh  nH 
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(1+M)(1-K){2+L-K)    sinh  nH      (1+M)(1-K)  ,       , 

: — ^;^T^2 hh 2rM+Kr  ^°^^  "^ 


+  -r—  sinh  nl 


T.(]±ria  3.•n^2   nH 
M+K      ^■'"^^       2 


Substituting   tlie    values    of  K  and  M : 


H 
f2 
T-i-  R   dn    =   A^    sinh_nH  _  ^ 

[j.nH        fj  '  n  ■(  n        nH  n  n 


*  ^(A-f2^)(l-e)l 
3A+2ti  J 

+   p2/2(AH-tx)(l-e)  L   ^  2(\+2n)(l-e)1si 


sinh  nH 

nH 

nh  nH 


-  cosh  n 


3A+21X 


1| 


,.   ,2    sinhnH_^   ^^       [i  ^  2^ +2|x)(l^)l si 
n        nH  n  n  U  3'^+2;j,        J 

+   d2/2U+u)(  l-_ei^  ^  2i|.2^1^n 


■\ 


nh  nH 
nH 

sinh  nH 
nH 


-  cosh  nH 


2(A+uO  1-e         ^      rr^mi   .    ,       „      [U+2u)  +  2\e]      .    ,2   n 
— ^TTTo 'Cosh  nH+-|— sinh  nH  -       -,-.  r: -^  sinh     -»- 


.'» Ozifiir. 


H 


[  J> 


n  '.1 


:  pt  Lt  X '{  :9.cr  V 


Aa***-  «.4  •» 


[p 


i        ^' 


Hr:    -1 


w 
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Examination  shows  that  the  terms  in  C   and  D   are  identical 

2  nH 
ith  those  in  ^      and  3  except  for  the  last,  in  which  sinh  — 

But 


^      J      ^  V.2   nH 

appears    instead   or   cosh     -y- 


v,2  nH        T     ^   „,-^u2   nH  .    „.„v,2   nH 
oosh     -5-  =   1   •*•   sinh     -^  >  smh     -5- 

Hence   any  pair   of   values    for  A      and   B     which  make   terms    in  them 
positive-definite  will   certainly  make   terms    in  C^  and  D^  posi- 
tive-definite   and    the    former  thus   give   the   stronger   criterion 
for  stability.      The  condition  for  stability   is,    therefore: 

.     sinh  nH  f2CK+ix){l-e)U   ,  2(  X-f2ix)  ( 1-e  Hsinh  nH     2(\-Hx)(l-e)         ,       „ 

nH      .    V      rr        [{X+2[i)+2Xe]  ,2  nH 

+  -^  smh  nH  -  ^^^^^ cosh     -^ 


1  + 


for  all  n.   This  simplifies  to: 


2{\+2\x){l-e) 
3A  +  2ti 


sinh  nH 


-,.2 


nH 


-  cosh  nH 


I 


[ 


f  ^    o„  Nslnh  nH  -, 
(1.2e)— ^^^ 1 


'I 


■  .  2(X-t-2^)(l-9) 


A  +  2|j, 


slnl'i  nH 
nH 


+  1   >  0 


which,  since  the  right  hand  factor  is  positive,  reduces  to: 


/,    o^^  sinh  nH  n  ^  ^ 
(1  -  2e;  -— -g 1  >  0 


or 


e  < 


nH  cosh  nH 


which  is  most  severe  for  minimiom  n,  i.e.  for  n  =  1,  giving  the 
limit  of  stability  or  critical  strain  as 


or 


1  -  H  cscb  H 


<  p    for  any  H,  no  matter  how  large 


H  ? 

-  j2        for  F.  «  1 


i    (CtSi;!    ;i>n;j 


:fj-fi 


w 


^Irt 


nnx 


'-  n  n 


0     <  !       :■ 


or)    y5.-i;i-. :    . 


-■«-  c. 


^■^    <  ^    — 


,;>    ,i    =   ii  ■10": 


S  '>  I     ;n  f .  '."  0  " 


■i"!;''i    t:,    yJ,  .;. '  d?.•■■:• 


£'^. 'TJ'?.l    '.■.'"'      ^" ;  J  r 
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2 
Furthermore,  e   «  1  implies  H  «  1  sc  that  the  last  approxi- 
mation applies  for  small  strain  as  well  as  small  thickness  com- 
pared to  length  of  colijmn. 

The  critical  load  is  the  load  producing  the  critical 
strain  and  is 

cr      X+ji.      cr     cr 
2    2  2 
■^^    12  l"^ 

for  small  h/L  or  small  strain  in  agreement  with  the  usual  Euler 
load  for  the  assumed  end  conditions. 

5.   The  Hollow  Cylinder. 

a.   Symmetric  Conditions. 

In  order  to  keep  the  problem  two-dimensional,  the  cylinder 
will  be  assumed  to  be  so  loaded  at  its  ends  as  to  prevent  any 
displacem-ents  parallel  to  its  axis.   Using  polar  coordinates, 
the  principal  strains  in  any  cross-section  are  given  by  equa- 
tions (6)  and  (7)  which,  substituted  into  equation  (3)  gives  the 
strain  energy  density  as: 

(8)   2w  =  ( X+n  )i/(iM^ +49  +  ^49)+  r^f*   [(pt7)ip-i  r^f  -  2 '  ^ 

with  the  total  strain  energy 


(9)  W  = 


..a  ,-.2% 
w  dv  =  I   1    wp  dpdQ    per  unit  length 
3b  Ja 


::i.,- 


".    ftriv    -liiri:;   o^   I   >>    "H   as^fqirs-t 


I  nm.i.' :  ■- 


.f-E:>i.t?To  oi'ii    -  .i.t  r^V'boxq   '.-^O'i; 


c    .        ^ 


.3 .  i  Pi:  e.:.;ir.~o%  '• 


:'0--i: 


-capo    YO   nsv.L:    ■^•■•le    ^  o  '  "r^>*^-f-ao-;  •■   yt.i,^   ;:";    ::i.."\"r:'s   I- 


/  - 


•Q'    o 


1  + 


iW'-;?] 


1  ,  I 


n:  ■*'■,■.  cj 


jt.n'!  •-■.■:>'\ 


13 
If  deformation  is  symmetric,  we  have  4  =  Tq  =  0  so  that 

(10)  Q  =  2Wp  =  (A+tx)p(^+r')^  +  ^p(|  -  r')^ 

2  P 

=  (X  +  2\i)^   +  2\rr'  +  (X+2ti)pr' 

where  primes  denote  differentiation  with  respect  to  p. 

Applying  the  usual  Euler  criterion,  W  is  a  minimum  for 

5r  ~  ^3r' ^ 

or 

2 

p  r"  +  pr'  -  r  =  0 

whose  general  solution  is 

(11)  r  =  Ap  +  I 

The  change    in  volume    at   the    outer  surface    is,    per  unit   length, 

Ti(a  +  Aa  +  -)      -  7ta^ 
a' 

making  the  change  in  potential  energy,  or  work  done  against  the 

pressure  p: 

V  =   7ip(Aa   +  -)(2a   +  Aa  +  -) 
'^  a '  a' 

Also,  substituting  (ll)  into  (10)  and  (9),  the  total  strain 

energy  is 

W  =  27r[(A+ti)A-(a^  -  b^ )  +  [iB^i^   -   —)] 

b"^   a 

so   that  the   total   energy   change    is 

T  =   W   +  W    =    7r(2(X+ti)A^(a^-b^)+2|xB^(-^-^) 
I  b        a 

+   p[2a(Aa  +  ^)    +    (Aa  +   |)^]] 


.J£fi''     02 


-       .   T       =■      .:•) 


I    • 


\    'i    -   --'^x^   t 


'.TC,(    ,S-«-/.  f     ■«-     -^.:l/S, 


TCC 


^c    -     '-r.:    -<    "-".^ 


noxrn  ^  -  3 


.riixiir.el    Tin:? 


-Ji'O    ►^^'■.■:     "' ?^    '^!:::\- Lr 


Qdn-   rt^rtl  A-^p    frtJiO 


'  ;:-.^-7  0-T    /: 


>r.Tr 


:i.CS'T^3    i^'^'O-f     ;:':io 


^XI) 


■"d 


(- ■:- 


i:  ^ 


>  :  .1- 


; ;:  >  , ::.    ■* 


lu 


This    Is    a  minimum  when 


i|.(X+^)(a^    -   ta^)A   +   p[2a^  +   2a(Aa  +  |)]    =   0 

and 

i|^(4   -  ^^^   "■  Pt2   +  |(Aa  +  |)]   =   0 
b         a 


t 


from  which 

(12)  B  =   iiL±.id  Ab^ 

and  v,2 

2(X+tx)(l    -  -2^) 


A 


(13)  P  =   -  ^ 


1    ^   A    -.  iAJ--El  .b^  A 


For   small   displacements ,    the   latter   approximates 

,2 
p  -  -2(x+^)(l    -  ■^■)A 

a 

in  accord  with  linear  theory.   In  (13),  however,  it  is  seen 
that,  even  for  infinite  displacement,  P  remains  boimded.   This 
probably  has  more  meaning  if  b  and  a  are  Interchanged  so  that  p 
is  an  Internal  pressure,  rather  than  external.   Then  we  hav?^  th: 
condition  that  the  pressure 

2 
2(X+|i)(-\  -  1) 

P  =  H — 

LL    2 

causes  infinite  expansion.   If  the  cylinder  is  thin  compared  to 
its  mean  radius,  the  maximum  pressure  approximates 


U^i(A  -^  ^x)  /a  -  h. 
P^    A  +  2tx   ^   b   ^ 


i  ;■•:■    }    I'.i  )£h 


^     ■   c; 


/) 


\-'    I  I.;    +    A  1 


+    .t 


rJpjviX'^'::ca>=    aert.ffii    i^rfj 


I)(.:;+X  C 


or 


:f  !:     ,-cyvet.-ori    ,  i  ::.i;    c^ 


f:''i^-i;:.&"i:f   '.ni^ 


IS. 


i.'  ;? 


t  ■  +   X 


■_-l,-J 
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b.   Total  Energy  in  Perturbed  State. 

We  now  assiime  small  changes  in  displacement  from  that 
occurring  ujider  symmetric  conditions  and  thus  take 

r  =  Ap  +  -  +  P 
P 

where  P  is  a  function  of  p  and  0,  and  permit  c|)  to  differ  si^^^^-t-- 
ly  from  zero.   Then  (10)  takes  a  more  complex  form  in  which,  if 
powers  above  the  second  in  P,  cjij  and  their  derivatives  are 
dropped, 

2 

0^  P    •  P    ^ 

+  [,.p[(X+^)A  -i^]P  +  (\  +  2^)p(^'-2+  P^)  +  2XPP 
P    ^  P     '         ^ 

+  2\p(l  +  A-.^)P^40  +  2[(X-F2ti)  +  (3X+i^a)A+^-^i3,]r,; 
p"   •  p 

p  p 

p   ^ 

The  work  done  agc-.inst  the  external  pressure  also  assumes  a  more 
complex  form: 


v  =  p  .[['■'  (£-1-1.)!  ,(e  .4)1         .  ^^^ 


VJo 

2%    (a   +  Aa   +  -  +   P)^  I  , 
^ i(Q  +  4)  L_      -   na^ 

27:  I  ^ 

P  n         F(a  +  Aa   +  I  +  |)d(9   +  4)  jp^^   +   7i(2a  +  Aa+ |)  ( Aa  +  |)i 


"■^  f 


^'■f.' 


:^£:;iJ    (r-"."Ll    Cf. 


Z,;.!:^.--       ^f,!' 


rr^      i_     ^"     -1.      ■    I*" 


J.    .x.-.^'>'1•  i-;^. 


1  ■  '•) 


Lii    -evfrh^-. 


' '.  ■") 


,^    £•' 


::,■:     ,h     ,"■     ns 


'   n©.ii.i 


•-  ..u;*-:  ]q.:^ 


■'!    ■;.'*■ 


;       -I  *  ~ . 

^*  '     ; 

1 


'.  1 


•"V    -    •■>      fi--       *:,:;.  t  i !  '■  ■;  n  '-'•      '!>  f ,  ; 


,-^:'.-r<..    i-\0' 


/• ; 


-If:     .<       ,:   V    ,  -(         1 
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But,    to  first   order,    for   p 


d(Q   +  4)    =    (1   +   cl-old© 


Hence,    to  second   order   in  P,    4»    and  derivatives 


V  =   Tipa^(A   +  -^)(2   +  A   +  ^)    +  pa      'f(1   +  A   +  ^)dQ 


Jo 


+  pa  P[|  +   a(l   -I-  A   +  ^^-)4^]d9 

Terms   in   Q,  and   V  which  are   independent   of   F  and   4   corres- 
pond  to   the    syrriLtietric   case    and    the    total    energy  correGponding 
has   already  been  minimized  by   our  choice   of  A   and  B   in   accord 
with   (12)'  and    (13)'      ^^'^   ^^s   n0T^/   turn  cur   attention   to   ter_ay    in 
Q  of   first   order   in  P  and   4    and    their  derivatives „      Sines 


2% 


,  d©  -  4 


271 


For   continuity,    the    term 


i^p[(X   +   tOA   +   i^](l   +  ^    +  ^)4 


2'^0 


does   not   contribute    to   the    strain  eaergy.      Also 


=    2p!:(?.    +    iJ,)A    -  ^i]F 


,2 
2(A  f  n)Aa(l  -■^)F 
a 


a 

b 

p=a 


in   view    of    (12) 


-=    -paF(l   +  A   +  ~)         in  view    of    (13) 


->  [,^:r 


.u'rv- ..'I'^s-L  •  r^-rn.    ,  i>    ^J^ 


;::.ivU 


'.     v    i.)( 


'■■^  •"  i      !  ■  '  .  — : 


re-icoo   lii   bfi' 


!.  -Nl  ■ ' 


.1   -^:-"T?>. 


■  I      C    .1        ; 


I :  .V.V       (^     ,"1  "JV: 


I   '-      .  l' 


••  •■■  S'la 


O    ^    r  1 


-1    ■■; » 


■x{v 


!    '  -.-. 


,..."'.'.   _     .'<  i 


1    /.  .■  :> 
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so  that  these  terms  in  W  cancel  the  first  order  term  in  V  leaving 
only  quadratic  terms  to  be  considered  in  the  total  energy.   It 
is  convenient,  at  this  point,  to  introduce  the  general  Fourier 
series  expansions: 

(11;)         F  =  P^  +  Y~  P  cos  nO  -  V~  f   sin  nO 
^'  o   >^—  n  ^ —  n 

and 

5~  £   cos  nO  +  '5~  \i/   sin  nO 

(15)  4  =  f,  -^  —-^ 3^:^S 

o  p(l   +  A    +  5^) 

P 

where   P   ,    F    ,    f    .    t"    ,    ^    ,    vj/      are    functions    of  o   ard    independent 
of   Q.      Also,    let 

(16)  K  =   ^i-i^S^  =  ^_  (1   ,  „   .  5j,    >  0 

if  the  center  is  not  to  close  solid  under  symmetric  loading. 
Substitutiiig  in  W  and  V  and  noting  orthogonality  conditions  for 
the  Fourier  terms,  ,;e  find,  using  primes  to  denote  differentia- 
tion with  respect  to  p. 


where 


/^w  =  f  P  R  dp 


R  =-■  R    +  \     r    +  3~"  s 

o        ^- —     n       ^^-     a 


with 

R^   =    2(A+2i^)-^  +   2(A+2[x)pP^^   +  l^XP   F^ 


1+K        PU+A+      2^^o 


VI 


■  S.   Y  n*  mni3:^  Ti'jb-i;; 


J I 


rf ... 


'xsi/ilfo'il   .LRi8ny:v;    orf-J-    cctrnotirtl    OJ    ("fiOq 


^.-i^   oz 


^a   r?:: 


c 


C'  C-'      '"'   ^  ? 


U5 


J.  /. 


^  -  <b 


<.;'; 


■i;..'  '  -  V-J  :^.    '- 


\j 


0  ;:t- t:''iff"^a    •^■J.:..;,.,     L,'!_.';.    -• ' 


:v::a. 


i  i-: 


'  ''i     'M  ■'.  ij 


■v.-^  ^ 


ii    -I 


i  /  -Q 
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R„  =    (X+2ti) 


n 


.p         ^  n 


2   2 

n  \j/  2nP  i)( 


n 


n^n 


M-M-ni^n^n  1    +  K 


2n  F 

n 


P       J 


+   2\F  P'    +   2\nP'\J/ 
n  n  n^n 


L     P 


2nP  \1(  p 

+  2nF  \|;'    +  S_il  +   pvj, '  "^ 


2(1  +  A   -  ^) 
P 


1   +  A 


B 


n  n 


(1    +  A    -   ^)2 
U"   A   +   ^?   "^^ 


S      =    R    (f    ,£    ) 


and 


Av  =   p  5  [^P"^   +  ^  (F^   +   2nF  i}/    )    +  y"   v^  ^    ~^^  €    )  ] 

2 
Sines  the  term  in  f '   Is  positive,  1a'  Is  minimized  for 

f;'  =  0  or  ^  =  0  if  simple  rotation  is  ignored.   The  tormg  in 
^o         ^o  '  ^' 

F  and  F''  are  independent  of  other  functions  and  correspond  to 
o      o         ^ 

symmetric  conditions  previously  treated.   It  is  easy  to  verify 
that  F  =  0  gives  miniinum  total  energy  as  expected.   Since  ter^ns 
in  f   and  §   are  identical  to  those  in  F  and  \1;  ,  we  may  consi- 
der them  as  mereljr  corresponding  to  a  rotation  and  limit  consi- 
deration to  the  latter  functions  only.   Vfe  note  that  functions 
of  a  particular  subscript  appear  in  separate  terms  from  those 
of  other  subscripts  so  that  the  total  energy  consists  of  inde- 
pendent contributions  by  terras  involving  different  subscripts, 
with  each  given  by 


n   2 


R  dp  +  p(F  +  2nP  \1;  ) 
n   "^   ^   n      n^n 


so  that  minimum  energy  corresponds  to 


5R. 


3R 


-2  =  (—2) 
n      n 


r. 


^•^1 


'^^^ 


[■      s,^  _ 


1  ,.':'    .:' 


(5^ 


(    iv  ■J^.-vi    ^ 


n 


•I' 


J.'^-*v; 


n- 


Pn^'.'   ^^^^    '^  n' 


n 
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and 

dR  aR      , 

which  leads   to   the   two   differential   equations : 

(l+K)(p^P"   +   pP'    -   P   )    -  Kn^P     =    -npij;'    +    (1  +  2K)n\l/ 
^     n        ^  n         n  n  ^^n  '    ^n 

and 

K(p^vl/"   +   pil;'    -   \^)    -    (1   +   r^n^il/^  =  npP'    +   (1   +  2K)nP^ 
n  n  n  n  n  n 

which  are  readily  solved  to  give  tbe  general  solutions: 


n    n^      „n-l    n*"^      ^n+1 
P  P 

,   -   r  o^+l  ^  (n-2)  -  2K  ^n    (n+2)  +  2K  j,   n+1    ^n 
^n     n         n  +  2K    n-1     n  -  2i:<.    n'^       n+1 

P  P 

except  for  n  =  1,  for  which  we  get 

P^  =  C^  +  D^  lop  p  +  E^p^  +  -^ 

*1  "  -^1  -  rr2-K  -  °1  l°g  P  -  1  -  2K  ^iP  ""  -2 

p 

Using  the  method  given  in  Appendix  A  and  noting  that 

.2(3_i JiJK     ^   g,^  at    p  =   b 

1   .  A   .  t^ 

P 

=   2\.i-p        at    p  =   a 
we    find   that 


:  Etoi  t£j 


,■■1  ^   ( 


a  ■ 


"U 

.T       ^ 


jn  =    _^^ir\i{yi    i-   i) 


■:;c   -Lf^'Torro- 


.  \ 


IX      -'     :^. 


Q      mOL       r'^' 


—  "*  J- 


tc::c»   ;?nl;f..i;-'    I:;-;    .•    ,.  "  :.ru '-rq-i   i\x    ■.:-,^^:  '^%  jjodi  z^f^ 


f'-jj-. 
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2T     /  ,r 

which  becomes,  upon  substituting  the  above  values  for  P  and  \|;  , 
(17)   !!n  =  2(n-l)c2[2^(a'^^-2  -  b^^"^)  -  pa^'^'^] 

2(n-l)D  (  o       ^       o-^    o        o-n   '^ 

+  ^-n      [2i(l+2K)[2(X+ti)K-^]  +  ix(h'^-l)J(a^^-^-b"^^-^) 

+  p|(n2-l)  -  (l+2K)2jb2^-2] 

^^"""^^^l  [2{(l+2K)[2(A+n)K-ti]  +  ^i(n2-l)|(a2'^''2  -  b^^^''^) 


+ 


(n-2K)2 


-  p|(n^-l)  -  (l->-2K)2|a2^-'2j 


2 

2(n+l)H         2n+2   ,,2n+2.  .  „^2n+2. 


2n+2,  2n+2 
a    b 


[;(n-l)(l  +  2K)  ^  ^  ^    kiIlillLl  +  2K)  g  .. 

+  ^(^^  -  1)  C  E  [2^(a2^  -  b^^'^)  -  pa^'^] 
n  -  2K    n  n  ^  ,        t-        ^ 


a  b 

p 
+  Q(^--l^  D  E  [(2tx(l+2K)  -  2(A+2u)K^'(a^  -  b^J  -  (l+2K)pa^] 

for  n  y  1,  while 


[{.^'^ 


It  Xj      n 


'i  1    -1 


K^u    .    ... 


.>.-    4      V, 


::  ~  rl  -3 


^  _. 


M    •'    i   (■! 


";  -J    .  '  .■> 


•        J. 


-.f:i  ^ 


I  ^ 


r        'I 

■-;   '    u;  ■-  "v-i; 


-•  <-  X 


•  -I-    i 


2,,.c 


q 


rr  f 


I  - 


t"/iqv:>ilii-i  /     - 


!:■  ^i  ;.Ui':    ;  .  i.';'. 


oinrff.'    ,f 
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2 


2T  D 

(18)      — i  = 5-— «   [i^(X   +  2^)K(l   +  2K)    log  f  +  p] 

^  (1+2K)^  ^ 


I|.E^(1+2K)       .  ,       .       ,  . 

— i 5    [U(>^+IJ-)K   -   2^L^(a^t■-b'+)    +    (l+2K)pa^] 


+ 


(1    -   2K)" 


+  Vt   [2tx(a^   -   b^)    +  pb^] 


a-'-b^ 
+  — i-i-  [(l|(X+ia)K   -   2a](a^-b^)    +   (l+2K)pa^] 

^^1^1  ro    .    2        ^2>    ^     ^2.        8(H-2K)    ^  „ 

(l+2K)a^b^  i    -   ^1^.        11 

and   is    independent  of  C,    as   might   be   expected   since    terms   in    C, 
correspond  merely   to    a   translation  without  strain. 

c.      Stability   of  Symmetric   State. 

For   the   symnetric   state   to   be   stable,    T     must   be   positive 
for  every  n.      For   convenience,    let  us   take 

(19)  k  =  ^  <  1 

a. 

1   +  A    +  ^ 

(20)  s   =    1    -  ^  = 1- 

a 

2(X+2ti)(l    +  A    +  5^) 

(  21 )  t   =    2(2^^2lx)K  ^ bl_ 

^  (A+3[x)    +   2\i.[A   +  —^) 

b 

all  of  which  are  positive  if  the  center  does  not  close  solid 
during  symmetric  deform.ation . 

Let  us  first  consider  T, ,  given  by  (l8).   We  substitute 
the  new  constants : 


re.. 


iq 


i)i  Ahli  V'~ 


5;:.  Mi. 


.'\x.-5  ~  i) 


s 


■'~)-'iS]     :;t^ 


■^  a"'  Si 


:>I) 


)   ■ 


x(!'+/;)4if] 


;,...^.. 


>•  -    ^-tfq    ^1    \- 


■;■    ■  i !  *-. 


.\-' 5'£  ^.  J  ■'.    .'co;::v;-v;    'lO  •  .  .r^.i.  .  rt  >  i    r    o-5 


■I;:;;-'    ':! 


[r-Vo      IC'      V 


k^  V    A.  u   J.  < 


•::■   -I 


■  i       W  .  t3  V  S3 


•»        f 


L  ; 


O 


H/; 


!-!.   'j 


t.i.ct;.     JOi\    Goofc    'l&^^o'::       .:.j     J. :     cv:  .t^  •••:■  V    OT^    .10 


.:t')r    ^■.''      .(bl)    "^cl  Li''vi.z    ,  .Z 


22 


X-,    = 


2H^ 

2(1    +   2K)ab   „ 
1    -   2K       "  n 


D, 


X-,   = 


3        1+2K 


and  note   that 


to   obtain 


logf  > 


a     -  b 

~2  2 

a"^   +  b^ 


-i  >  x?(k'^   -  sk^j    +  xf[(t-s)k"^    -    (t-l)k^] 


-f  x^Ltd   +  2T)^  ii^ 


-lil  ^   1    -   s] 


k       +  k 
+  2x^X2(1-3)    +  2x^x^(k"-'--sk)    +   2x2X^[(t-s)k 


-1 


(t-l)k] 


Now  the   quadratic    form 


y~  'S     a.  .X.  X  .      ; 


^ij   "   ^jl 


is   positive  definite    if 


Ai   =    a,,   >   0 

A2  = 

^12 

^12 

^22 

>   0 

'-^11 

^12 

0 

■  13 

A3  = 

^12 

^22 

^23 

^13 

^23 

"33 

>    0 


etc.   In  the  present  case, 


o .: 


•:S.-K 


-J      -r  li 


•>  \ 


.>^  V   ^.■■•'■ 


r  I 


\       .4  , 


r  ^  ^ 


■«    1 


-ji 


-n  ?• 


v<; 


(.  'j  ._•  - 


(2- 


■-■f.-,-'      •-'J 


n- 


i.  ^  .... 


"'  < 


.1,  L 


i  r' 


)    ^:-:'' 


L^.:i 


;     .1 


\ 


i 


•■T  ^ 


,<„      , 


:-'\ 
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At  =  k-2  .  sk^  >  1  -  sk^  =  i^^^ni  -  k-^)  ,  0 

^  1   +  A    +  ^ 

a 


A,  = 


-2        2 
k  ^-sk^  1    -   s 


(t-s)k"^-(t-l)k^ 


1-3 

=   [(t-s)k"^   -  s(t-l)k^][k"^  -  k^j 


but 
and 


-2  2 

k  ^    -   k      >  0 


(t-s)k"^   -   s(t-l)k^  = 


(1   +  A   +  -^)M 


(1   +  A   +  ^)[(A+3^)    +  2ti(A   +  ^)] 
a''  b^ 


viherQ 


M  =    (X+n.)(l   +  2A   +  ^)(k'2   -  k^)    -  2(X+2ij.)^„-(k"^   -   1) 


=  -^J^    [^t(k"2   -  k^)    +   2(X   +  2t.)A(l   -   k^)] 


[2(X   +    |i)K(l    -   k^)    +    ti(k'^    -   2    +   k^)]    >  0 


=  iiiil   [2(X  + 
Hence   ^.2   ^  °'      Finally 


A,  = 


-2  2 

k  ^  -  sk"^        1   -   s 


k'-'-    -   sk 


1    -   s 


(t-3)k"^  -  (t-l)k^        (t-s)k'-'-  -  (t-l)k 


k~^  -  sk  (t-s)k"^  -  (t-l)k  t(l+2K)^    ^^_^'^1+  i  . 


-1 

k"-'-+k 


(k-\kf 


-2 


-2 


sk"^  +    1 


k  "    -   sk"        1    -   s  k 

1-3  (t-s)k"^-(t-l)k^         l-s+(t-3)k"^  -  (t-l)k^ 

k"^-3k^+l_3       l-&»(t-s)k^-(t-l)k^        t(l+2K)^(k'^-k^)  +  (l-s)(k~VK)^ 


{->l 


-    I    <     ">iS    - 


'i((  r-^)-^~^(£^ji 


t   1 


[- 


^... 


>!   ![■   :l{f~i} 


m^-  ■<  A  -c 


^^->. 


i  «•"■ 


A  +    £) 


I- 


;,  •"   \ 


AS    4-    j 


i  ( 


DA'.-.  ^    /:;?   +    '^>'  -   ^'^fV-J   ■'^^ 


s- 


'■■''■       L  ', 


X  J  I 


r  '       ,        J-  '' 


,0   < 


■*  ,  r         .^ 


_    r 


.,  t ..  -t  \ 


■si 


(  V-rf 


:|,  -    ^~:i^e-J) 


>*-' 


\  •■-■ 


•■■  I  ?  ■■ :! ; 


L- 


i      I 


.;;<■-::}  --"'L^^-rO 


M)!c.  -i^  ■!• 


>:)-'Ui^'.(: 


J  '^Z  • 


21^ 


(k--^  -!■  k)' 


-2  2 

k  '  -  sk    1  -  s 


0 


1  -  s 

0 


(t-s)k"^  -  (t-l)k^   0 


^2' 


(k"^  +  k)2 


where 


L  =  t(k"^  -  k^)[(l  +  2K)^  -  1]  >  0 


so  that  /A^o  >  0  and  1\  is  positive  definite  and  is  a  minimum 
for  D^  =  E^  ^  H-j_  =  0. 

We  now  turn  to  T  ,  as  given  by  (17),  making  the  following 
substitutions,  as  well  as  those  given  by  (19),  (20),  and  (21): 


n   = 


n^  -1 
(1+2K)^' 


0 


2(n_+  Ijy^  ., 

il+l  a+1        "'n 

o"2~'  ,  2 
a    0 


^2 


n-1  n-1 
2(n  -  1)^/^  a  2   b  2   C 


n 


^3 


1  /y   H"*"!   i2"^i. 

2(1  +  2K)(n  +  1)-^^^  _  2   ^  2 

n  -  2K 

2(1  +  2K)(n  -  1) 

n-1   n_-l      '  a 

(n+2K)a  ^     b  '^ 


1/2 


a 


D 


E 


n 


to  obtain 


u 


p    ~ 


"v(I-^  ^ 


5i-   , 


i;\2-:t) 


,"»     -      1 


'(-'■■    ^    f  ;  :v*'>(   ~   -^'^J);*    ^ 


nriirh Jh!  b    ai   fc;"^^^   sr^Znll 3xj  svJjir...^. 


'vOlIqI   sri:*   i,-i.oI 


:.  ■'-    . .  :.'   o:l"   :T~;f<7 

•       ,:&     S..i      W^     ^.B      , 


.1  ?-.o         i   Ci 


n- 
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7111  1  ■  Z 

+    x|[(t     -     3  -     S,;^)k-^-^     -     (t     -     1    +    //2)k'^"'^] 

+  x?[(t    -   1  +    ^^)k"''"^^   -    (t    -   s   +   s/)k"'^] 

-  2x-,x-,(l    -  s)    -   2x2X>  ( 1   -   s)  •  ■ 

+   2x^x,  r](k'^   -   sk^)    +   2x2X^>?(sk""   -   k") 

-  2x^xj^7/[(t  -  2s)k"^   -    (t    -   2)k] 

Now 

,  -(n+1)  ,  n+1   ^  ,  -2  1  2   ^  ^,  ,     „ 

k  -  sk  >  k        -   sk     >  0  as   before. 


Further,  an  obvious  criterion  for  stability  is 

sk"^^'^-'-  -  k^"-"-  >   0       for  all  n  >  1 
so  that 

s  >  k^ 


Assume  this  to  hold  and  make  the  substitutions 


a3_  =  k-^-1  -  sk^^^-  >  0 

ag  =  sk-^-'l  -  k^-1  >  0 

b^  =  k"''^  -  3k^ 

bg  =  sk-"  -  k" 

b,  =  (t  -  23)k"-'-  -  (t  -  2)k 

1  -n-1   ,  n+1 

c,  =  sk     -  k 

„      _  ,  -nH     ,  n-1 

Cp  =  k     -  sk 

d,  =  (t  -  slk-""-^  -  (t  -  Dk""^^ 


d2  =  (t  -  Dk"^^""^  -  (t  -  s)k''-^ 


1+n- 


Li-n  ,  l-a- 


Jla)^x   +    0        'Aq   - 


[^^•«.u^  .  ..  -  yy 


5i  ^    \^  f;    I-    a    -     ■'  /    -    ■  ••.    '.;  V I     ;'.    +    X 


(E.-.i;_^,x...i^.-   (5  -■ 


!•" 


i'jixS  -    ^)   -   ■'"■••isS 


xS 


. otoleo 


hv-  .,     ,   I-in 


r.  vcff  J.  •  cictft  -Tc''    ::fo  !:".   '.r  ir'  Sj/o.' "vdo 


I  <  n   .£/ 


i-.-t  ■        X-s-n-  , 


:a">  -'c.-i/.y  a:. 


^?>j;r-    _:c\'-.    II-;)--.;  .'l'-:^ 


r, 


-■•■■f'-.'. 


f  .    „ 


r:      .    f; 


■1  ] 


■f; 


>  t; 


-r-;: 


,!.-'  a- 


^^:;('; 
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Then 


2    .  b^ 


2   2      r.  \2   ,b2 

.   2r^  .2  ,  ,     (1-s)      Z      1-1 

+  X^[C2>.'   +  d2  -  — ^  -  >/   — ] 

-  Zvix^x^ih^    -    (1  -  s)^  -  (1  -  s)^-] 


and  is  positive  definite  only  if 


P  b^        b,  p 

,2(^3  -  (l-s)-A.  (l-s)~l2 


,U,-il^'i-,^(|-=,)!U,-<i^-.A!i-c,„ 


or 


)7^j[a^b^  -  2(l-s)b^][a2b^  -  2(l-s)b2)  '■"  (l-s)^(a^c-^  h  a2C2  -  2b^b2) 

<  [a^d^  -  (l-s)^][a2d2-  { l~s  )2]  + '^^(b^  -  a^C2)(b|-  a2C^) 

or,    inserting   the   values    of   the   a's,    b's,    c's,    and  d's   and 
factoring; 

[.^2(^-1  _  ^)2  _  (^.2n^  ^2n_  ^-2  _  ^,2 ,  ^  |,^2^2(  ^-1  _  ^^2 

-  [s(t-l)(t-s){k"^-''+k^^)  -  (t-s)2k-2-s^(t-l)V]|    >  0 

Now  k:!L-   k!!  .   k-n+l   ,    ^-n^3   ^...^  l,n-3   ^  j^n-l   ^  ^ 

k"-^-k 


nsfir 


^  .i-^ 


(i-.n 


r>     -f 


t;. 


vf        •»■ 


i"  -^  -^  ^'  ..'.^       :^-.V     X  h      4-      *"' 


^^ 


.r  s 


r- 


+    -^.'..o 


f^....- 


.     -    i)    -     , 


lij    '.'.■;.-    &S£i.r:  i'b  e--:       .  .  ^'i 


o  ?•    f> 


.L..(g.X) 


7^ 


L( 


■\  f. 


\     .     -" 


.0  ;    > 


:.••-'■}  ■!■  v-'^'.v^- 


's  J. 


r.        +         J  <-         -J        ^ 


\j;^C,3  -  c;-')^C 


;  -»  i"(;vi)  -,  Lvc;'''(?-X)  >   . 


bn.R   3'b   b 


0    .c'A    , -=; '  f  ■    sr'j    ^^    kw.).'..'^' 


■>;  +■ 


^;        *    4-  "^  *^'    ; 


-).    I  -:;-:(;-  ^ 


^f  •^ 
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Hence 
2,,  -1   ,  ,2      ,,  -2n  .  ,  2n  ,  -2   ,2 


^ 


(k"-"  -  k)^  -  (k"^''+  k^"  -  k'"  -  k") 


=  /;2(k-l.k)2,  (k-^-k^)2-.(k-l-k)2 
<n^(k-^-k)^  -  (k"''-k'')2  <  0 


Hence  stability  requires  the  second  factor,  in  the  above,  to  be 
negative,  or,  rearranging  it, 

(n^-  l)(k'^-  k)^+  [(^-  1)(1+  2K)k'^-  (t-  l)(l+2K)k]^ 

<  (t  -  1)(^- -  1)(1+  2K)^(k"^- k"")^ 

Substitution  of  values  shows  that  •  _  '  , 

(--  1)(1  +  2K)k"-^  -  (t  -  l)(l  +  2K)k  =  k"-^  -  k 

making  the  criterion        ^  . 

n^(k"^-k)^  1  (t  -  1)(|--  l)(i+  2K)2(k"''-  k"")^ 

Since    (k"'^-k")/n  increases    with  n,    requireraents    are    strongest   for 
n  =   2,    for  which 

[^.(k"^  -  k)^   <   (t  -  1)(-  -  1)(1  +  2K)^(k"^  -k^)^ 

or 

[1   +   2U+2ial  ^^j^      ^   2{M^^  ^^2]    ^    (2 j2 

^  ^  -     k-1    +   k 

or  P 

^  =   [^(A±2ixj2^2^2  ^   2(^-^)(l   +   k2)A   +    (l_lJS-)2   >    q 
IJ-  I-''  1    +   y^^ 

which  is    obviously  satis  I'ied    for  A   >   0    (    p  <     0)    and  obviously 
not    satisfied    for 

2lli2iil  A   =    -1 


oneli 


s 


■>i  -  '^'•'-l  f 


!,) 


J  (•  n. 


■\i 


V 


"i  r  n  > 


;^ 


(.. 


-•>:{}i^  f  D^.i  -  cf)  -  '">i(  'iS  •!  .r  >  u: 


'(3^ 


•"--  -f  J 


f  ;  I    ;■ 


,-y  1  i  ••  .-f  y 


I'^^jul    nicorir    ct^i.'/nv 


A  -.   ^-M 


^. ',,  ^1^    I     '.   I  I  J.   •-   J 


-"MU'S  ■;•  LWI  -^) 


s    • 


-•xoU    f%.^'-'\^xo'i^''-    '?"' : 


/\--    -  ■-/   V 


-  «! 


•/  -■.   - 


J;        r  ia 


t.U- 


■':h    ■••    .f 


■'i 


r-    .-.    A 
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corresponding  to  closure  of  the  center  to  half  its  original 

size,  since 

b(l  -f  A  +  ^..)  =  b[l  +  ^^  *   ^^"^    A] 
b^  ^ 

Thus  instability  always  occurs  before  compression  is  this  large 
for  any  value  of  k.   Furthermore,  since  s  increases  monotoni- 
cally  with  A,  we  have,  below  this  limit, 

s  .  k^  =  tLJ:_lAl2i^lA   /_  ^2 

\i   +    [^    +    (\+ij,)k  jA 

>.  +  2ll  ,2 

> c ^  _  1^ 

2\   +  3[i   -    {X   +   |x)k 

(X  +  2^)    +  (X  +  ^)(1  -  k^) 

so  that  instability  due  to  negative  ^  always  occurs  at  a  lower 

2 
deformation  than  that  corresponding  to  s  <  k  ,  and  the  limit  of 

stability  is  given  by  the  less  negative  root  of  the  equation 

A,  -   0,  namely  by 


^  1  +  k"^ 

For  1-k  «   1,    this    reduces    to 

2(X   +  2.0    ^  ^        (1    -   k)^ 
II  2 

and,    from   {13) , 

M,(X   +   |i)     ,,         ,  x3 

in  agreement  with  other  methods.   V/e  may  also  note,  that  |A  |  is 
the  average  between  radial  and  tangential  strains  and  is,  tnore- 
fore,  small  when  strains  are  small.   Examination  of  ^  indicate.?? 


r6nJ:g.i-T:o   e,-.^l  1j.nd  oi   -t'-l. 


1^1.:L±L  -,    r-,^  ^ 


!       A 


I)J 


■. :  X    zli!^    ai   no '•Ese-T.^i'-oo    a'T.o'ir^o'    •^.tA'oco    svjfc-;.:.*'^'' 
"/:icc)onoxn  £:aBj5e^xr.ri !;   ;^    of^'n    .O'l^^rm  C"lJ''-?f.'^,      «yi 


'•,.>,'       f   r   -1  +      :.r 
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.w 
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itr  _."!  _:,  _ < 
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n. ;  n.j 


f     .-I 


/;tl 


07    li  "^-^i  b's'i    i.'  :!,t    ,  .■    ' 


ei    I  ^)     J/^fi.r    ,  sd-O;;:    cp  ■;; 


-  C-.     p  n  t- 
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X! 


^  \  f  i  }   n^o-r':.    ^-^u-'r 


I ..  <■'. 
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-that  Instability  can  occur  lor  small  |a|  only  if  1-k  «  1, 
hence  the  classical  formula  is  vn.lid  for  small  strains. 

6.   The  Hollow  Sphere. 

a.   Conditions  with  Pull  Symmetry. 

For  axi ally -symmetric  deformations,  eouations  (6)  and  (7) 
are  valid  for  the  principal  strains  in  a  plane  passing  through 
the  axis  and  the  third  principal  direction  is  the  tangent  to  a 
line  of  latitude.   The  change  in  length  of  a  circle  (see  Pig.  5) 
coaxial  with  the  assumed  symmetry  is  in  the  ratio 


dS 


3  _  (p  +   r )  siu  (0  +  d)) 


dS         p  sin  9 

so   that    the  third  principal   strain  is 

dS- 
e,   =  Yc<      -   1   =    (1   +   --/(cos   (])   +   cot   0   sin   (|))    -   1 

and  the  strain  energy  density  (l)  becomes 

(22)      2w  =    (X+ii)[2  +  {^-+  cig+|-  4Q)+rp]^>  (  A+li)  [  (  P+r  )4p  -  |  r^]^ 

+    (X   +  2^)(1   +  ■g)^(co3   (j)   +   cot   Q   sin  c^)^ 

-   2(3?^  +   2txj(l    +  ■~)(cos   4   +   cot   Q   sec   4)    +  3{3'^   +  2[i) 

+  2[\(1   +   g)(cos    (1)    +   cot    9    sin   4)    -    (3^   +   2ii)  ] 

r 


't2   ^-    (-^  -  4e   +  I  4^)'-   rp?T7T7^   r)4p   -  I   r,]2 


rc 


s  li/cf!--:  ul 
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and  the  total  strain  energy  Is 


W   = 


f[ 


,a  ^.r. 


w   d.v   =    2t[   \      \       wp"    sin   9  dpdO 

JbJo 


For  fully   symmetric   conditions  j    4   =    r^  =   0  and,    using 
primes    to  denote   differentiation  with  respect   to   p, 

Q  =   2wp^  =  U(X  +   ^)r^   +  UAprr'    +    {X   +   2ix)p^r'^ 

The  strain  energy  is  a  minimum  if 


or 


c^r  "  V^r'^ 


p^r"  +  2pr'  -  2r  =  0 


whose  general  solution  is 


r  =  Ap  +  :.^ 


P 


which,    substituted   ?.n   the   preceding,    gives 

W  =   2-k[{3\   +  2m,)  (a^    -   b^)y\^   +   liii{-^   -  ^)B^] 

b^        a^ 

The  work  done  against  the  external  pressure  is  the  product  of 
pressure  and  change  of  volume  or 


V  =  !i  iipa-^  [(1  +  A  +  ^^)3  _  1] 
-^  a 

making  the  total  energy: 

T  =  2Tt((3X+  2ti)(a^-  b^)A^  M^tx(-~-i-)B^  +  -|  pa^  [  ( 1  +  A  +  ^)^  -  l] 

b^   a      ^  a^ 

This  is  a  minimum  when 
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and 


(3X   +  2p,)(a^   -   b^)A   +  pa^d    +  A   +  ^)^  =   0 

a 


^         ^-)B   +  p(l   +   A   +  ?-^)^  =   0 

a 


L 

^"^?  ■   a3 

from  which 

(23) 

B  = 

and 

(2i|) 

P   = 

b3 


(3X    +    2p.)A(l    -    -i^) 
a^ 

U    ^   A    +i3^±2Ml    .    b3         2 

Inspection  of    (2Li.)    shows    that   p   increases  with  A,    reaches 
a  maximum,    then  decreases.      This   has   more   meaning   if  b   and   a   are 
interchanged   to   obtain   the   case   of    internal   pressure: 

.3 


P  = 


CU   +   2ti)A(-^'^    -    1) 


[1    +   A    ^  il^iiO    .    4  A]' 


whose  maximum  is 


3 

(3^  +  2'^)(^  -   1) 
b- 


and   occurs  when 


I      f    ^"Li'ii        sl 

A[i   +   13X   +_2iil  a_2^    ^    ^ 

for  which  one  easily  verifies  that  the  inside  of  the  sphere  has 
expanded  to  twice  its  original  radius.  Hence  an  internal  pres- 
sure sufficient  to  expand  a  thin  sphere  to  twice  its  dimensions 
will   cause    infinite   expansion  and    in   fact,    less   pressure    is 
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required  to  maintain  such  expansion  the  further  one  goes  in  the 
process.   This  effect  can  be  noted  in  blowing  up  a  balloon. 

b.      Total  Energy  vjith  Axially-Symmetric   Perturbation. 

We  now   consider    sraall  displacements   from  the   fully- 
symmetric  deformation  by   allowing  non-zero   values    for  4   and 
taking 

r    =    Ap    +   :§_    _j.    p       -   ■■ 
P 

with  (|)   and   P  functions    of   p  and    0.      For  convenience,    instead   of 
cj)   itself,    we    take   the   new   variable: 

\|/   =    p(l    +  A    +  ^)4 
P 

Substituting  in  the  expression  (22)  for  strain  energy,  we  obtain 
to  second  order  P,  ^,    and  their  derivatives: 

2 
2wp^  =  3[(3^  +  2p,)A2p2  +  kiK   ^-] 

P  .  . 
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+   2XpF   \|(   cot    0   +    (X   +   2(j,)\1(q   +   2X\1/\1/q   cot   0 

2[(X    +   ix)    +    (5X    +  kv^)A   +    2(X   +   2ia)-j] 
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+  2(2u   -   3XA)         r 


2   +  2A 


B-  ^e^P^p 


(1    +  A    -  -^) 


1   +   A    + 
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,B 
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+    (\   +  2|i,)\[f^   cot^Q   - 


[(3^   +   2p.)A    +   2ii  ^] 
pf_    ,2 


1    +   A    + 


Taking   as   element    of   volume,    the   conical    shell  with  vertex 
angles    Q   and   C+dO,    we   find   the   work  done   against   external   pres- 
sure  as 

V  =  ^  pi       (a  +  Aa   +  ~  +   P)^   sin(Q+4)d(Q+cj>)  I    _      -  ^  pa^ 


a 


P=a        3 


But,    since    (t  =   0   at   9  =   0   or  it. 
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It 


J  0 


(a   +  Aa  +  ^)^   sin    (Q+fj))   d(0+c|)) 
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1^    (a   +  Aa   +  %)^ 


Also,    to   first  order, 

\|;    cos    Q   +   ij.'^    sin   0 
sin   (Q   +  4)   d{Q   +  4 )  L_      =    [sin  ©  + i ]dO 
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so   that,    noting   the    relation 


P(l   +   A   +  ^) 


p=a 


^TT 


P(\l/   cos    Q  +   \!ip,sin   0)d9   ==    Pv};    3  in    0 
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we   have,    to  second   order, 


V  =  ^     paMd   +  A  +  ^)^    -   1] 
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+   2Ttpa'^(l   +  A   +  ^y 

a^       J  0 


ir 
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P  sin  0  dO 


p=a 


+   2iipa(l   +  A    +  ^)   [       (P^   -   P  4,)    sin  0  dO    I    _. 


a^     J  0 


The    zero-order  terms    in  '/.'   and  V    correspond   to    the   fully- 
symmetric   case    already  considered.      We   now   turn   to    the   first- 
order  terms.      V/e   have 
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in  view  of  (23).   Comparing  with  (2)4.),  vje  thus  see  that  first- 
order  terms  in  W  and  V  cancel  as  expected,  leaving  only  the 
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quadratic   terms    in  P  and   i|f   contributing   to   change    in   total 
energy. 

Considerations    of   continuity   require    that   \|;  =   P^  =   0   at 
0  =   0   or  n.      If  we    replace   9  by   the  new    independent   variable 

.•:  X  =   cos    Q         . 

and  take 

P  =   P(p,x) 
and 

\}r  =   sin  0   G(p,x) 

these  conditions   are    satisfied.      The   previous   formulae,    retain- 
ing only   second-order  terms   and    inserting   the   value   of  p,    then 

become  J 
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1 
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Sf 
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the   latter  may  be   written   as 
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We  now  introduce  a  new  function 
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If  (27)  and  (28)  be  used  to  express  G  in  terms  of  F  and 
this  be  inserted  in  (25),  we  get  the  differential  equation  for  P 
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tions   to   fully-symmetric    (4=0   and  F  independent    of   0)    for 

which  the   unperturbed  deformation  minimizes    total    energy  so   that 
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Appendix  A   Is    reduced   to   a   single   Integral,    given   by 
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c.   Axlally-Symrietric  Buckling  Criteria. 

We  confine  attention  to  compression  under  external  pres- 
sure, that  is,  to  A  <  0. 

Examining  the  expression  for  the  n-th  component  of  total 
energy  (3^))  we  see  that  terms  involving  J  vanish  for  n  -  1  as 
expected  since  this  corresponds  merely  to  sliding  the  sphere 

along  the  axis  of  symmetry. 
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For  n  >  2,    the   coefficient   of   J     must   be   positive   if   the 
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Prom  (32),  we  see  that  the  upper  limit  is  reached  for  a  lower 
absolute  magnitude  of  A  than  needed  to  make  S  negative.   We  may 
thus  take  the  criterion 
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so  that   Instability   Is    alvjays    reached  before   the   external  pres- 
sure   is    sufficient   to    reduce    the    inside    radius    to   one    fourth  of 
its    original   value    (lower  bounds   raay  be   determined   by   later 
steps    in   the    analysis). 
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or 


( 1 +^S )  u  >  £i3ii>  i_2il3_L-lM 


which  is  a  transcendental  inequality  in  S  v?ith  terms  dependent 
upon  n  and  the  ratio  of  the  elastic  constants.   In  particular, 
if  [J,  »  \,    corresponding  to  zero  Poisson's  ration,  the  lowest 
limit  for  S  occurs  for  n  =  co  and  is  S  =  l/2 .   For  X  »  [i,    cor- 
responding to  incompressibility  under  hydrostatic  pressure,  or 
Poisson's  ratio  =  l/2,  the  limit  for  S  is  1,158  and  occurs  for 
n  =  3.   In  the  former  case,  we  find  that  A  =  -l/3  and  the  center 
of  the  sphere  cannot  close  beyond  half  its  original  radius  with- 
out instability.   In  the  latter  case,  A  -  -0.i|66  ^  and  instabili- 
ty results  if  the  inner  radius  is  reduced  below  0.651  of  its 
original  value.   We  thus  conclude  that,  in  all  cases,  instabili- 
ty will  result  before  external  pressure  is  sufficient  to  halve 
the  inner  radius  of  the  sphere. 

We  now  turn  our  attention  to  the  case  of  a  thin  sphere 
and  let 

(1|9)  y  =  1  -  k  «  1 

2 
We   will  find   that   3    is    of   the   order  y   and   n     of   the   order   1/y 

for   the   limit    of    stability.      Keeping   this    in  mind   to   insure 

retention  of  all   significant    terms,    we   find: 

^  ^^/v2n+3)[l   -    (n.l)y   t  LnllU^ll  /]    +   6S  \ 

^  ^-^r(2n-l)[l   ^  ny   +  nl£n-LJd  /  ]    .   6S '^ 
n2        ^"-^  V  ^  J 


:'^n 


p;:!'.r   'I.-:'    Ic'-i-O   '^ 


''.^!     Of--Oi!i:>0-.C    £J.     ^'.rrc.DX?!    'IwilnJ     a. 


-■     ,f^^-njs?    I  >  r\    w.i.    ,T;:.i:':-    ■■•■:.,rj_f,.;  co    tiT:    f-y-:' 


X):  .."■:  .'J,,  e      j: 


• .'  ;  .1.  f   !; 


91 ana e 


WOiI     c 


'>  >i   -    i 


^;ii-l   X. 


+  r< 


■ii 


'■     -     flit  •'■.-C:  , 


•    f-fsj^ 


■£    -«- 


51 


2k2    [^    ^=    1    -    k^   =    y(2    -    y) 


J  k  X      *^ 
(2n-l)   C   x^""-^  dx  ^  (2n-l)y[l   +  ny   +  ""^^^^   ''   ^^   y^l 


2k^  [        .^   '^^      .      ^y(2   -  y   -   2S) 
k  x^(l+Sx^) 


k  x2^^^+Sx^ ) 


(2n+3)k^"^   1      -gH+r;:^ J-   -^  {2n+3)7[l-S-{n+l)y  + 


^'   x2-2   ,,  ^^  (2^.1)    ri    .   3   -.  ny  +  ^2lU     2j 


^^ ^^  -^  (2n-l)y[l    -   3   +  ny  +  ^^^- 

1  4   Sx^  ^ 


(2n-l) 
so   that 

+  3n(X   +   2ja)S: 


+  J(n  +  1)(\  +  2n.)s| 

t  in+liiinll)   ,^2]    .  3(/   .  2„  .  2)(X  +  2^)s] 

1^  -    Pn2  '^  ^nrrifeTTT  [[(2n^l)X.2(n2.n+l),](l.ny 
+  ai^lii  y2]    -   3(./  +  l^n  +   1)(X   +  2n)s| 


^)y. 


«      f 


.  < 


im)    -    L\'i{lJ^i&)  ^ 


■f-    III  -v   J  iYd-n;-  ? 


oii. 


f  p/: 


S)t.'-' 


r'^ 


•'i:^J.U;!iLL  ^  v-;  1  f^o-a-x ]vi c-' -^ ) 


+  -n\   ■■■'    I   -    Cr.:(  f-.!S 


1  .- 


'->:e  i-  I 


^  if  . 


■:-S\^J^C\\'i 


;,  (_   £•!    i^S  }.!< 


i^;  I^- 


i^is-j '; 


r  \ 


c  ^  ;;S    -r    A)[  J     •<•    r:  / 


T^-lO  r^  !'•!{;     m!- 


1  _•> 


ia^' 


)?^(a^S   •»-   a)(s    -   nG   - 


■•?;■ 


vn  - r  U  Ji^  ( X  !  r  •^■^xi ;S+ A  ( I-  "^r  } 


52 
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or,    dropping  y   compared   to    1,    and   higher  povjers    of  n   , 
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and 

p  (critical)  =i^y^.^3|||i^(l-.3, 


which  agrees  with  results  obtained  by  other  methods. 


Appendix  A 

When  an  integral  is  minimized  by  the  application  of 
Calculus  of  Variations  to  functions  in  the  integrand,  the  value 
of  the  integral  for  the  minimizing  functions  can  be  obtained  in 
a  simpler  way  than  by  direct  substitution  and  integration,  par- 
ticularly when  the  integrand  is  a  homogeneous  quadratic  express- 
ion in  the  functions  and  their  derivatives.   Consider,  to  be 
specific,  a  double  integral  involving  two  functions  in  this  way. 
The  general  form  is : 
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Now   Q  is   a   homogeneous    quadratic   expression,    hence    it  may   be 
written  as 

2  dP  X   op  Y   c^F  c»G  X   o'G  y   3G 

Now,    integrating  by  parts: 


p     3Q_  dx  =   P  5^ 


d       Jd 


f°   P  -4   (!§'-)    dx 


Pjx    '?P. 


f^  c 


"  d 


G 


9Q 


%-  dx  =   G  ^ 


3Q 


X 


G  4    ($^)    dx 
d        ^'^     ^^x 


a 


^y    3P      '^J'^         ^    9F 


pa 


3  b 


"^       y 


>a 


r:)Q 


c?Q 


jb     -y        y  y 


^. 


a  .-'  a 

-  G  #-  (|^)   dy 

,  ,          3y  3G           •' 

b  J  b           -^  y 


so   that,    in  view   of   the   minimizing  equations, 
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strain  Alon^  Principal  Directions 
Cartesian  Coordinates 


x-dS  sin  a,y+dS  cos  a 
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strain  Alonp;  Principal  Directions 


Polar  Coordinates 


p+dS  cos  a,© — i  sin  a 
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Strained 
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Buckling  of  a  Column 
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The   liollow   Cylinder 


Initial  Condition 
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?he  Hollow  Sphere 


^ xis  of  Symmetry 


Initial   Condition 
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Symbols 


approximates  to  the  specified  order 

X,\i  Lame  elastic  constants 

E  Young's  Modulus 

a,b  Radii  of  cylinder  and  sphere  surfaces 

k  Ratio  b/a 

L  Length  of  column 

h  Thickness  of  column 

W  Total  strain  energy 

w  Strain  energy  density  per  unit  vol^ome 

V  Total  work  done  against  external  forces 

T  Total  change  in  energy  of  system 

P  Total  force 

p  Pressure  per  unit  area 

e  Unit  strain 

dS  Increment ial  length 

n,m  Integers 

x,y  Cartesian  coordinates 

p,0  Polar  coordinates 

u,v  Displacements  in  x  and  y  directions 

r,4  Displacements  in  p  and  0  directions 

a,  (3  Angles 

/\W  Change  in  W 
P,G,\l;,f,^   Functions  of  coordinates 

Q,r),x  Auxiliary  independent  variables 
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